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Solvent plays an important role in the relative motion of nanoscopic bodies, and the study 
of such phenomena can help elucidate the mechanism of hydrophobic assembly, as well 
as the influence of solvent-mediated effects on in vivo motion in crowded cellular environ- 
ments. Here we study important aspects of this problem within the framework of Brownian 
dynamics. We compute the free energy surface that the Brownian particles experience and 
their hydrodynamic interactions from molecular dynamics simulations in explicit solvent. 
We find that molecular scale effects dominate at short distances, thus giving rise to devia- 
tions from the predictions of continuum hydrodynamic theory. Drying phenomena, solvent 
layering, and fluctuations engender distinct signatures of the molecular scale. The rate 
of assembly in the diffusion-controlled limit is found to decrease from molecular scale hy- 
drodynamic interactions, in opposition to the free energy driving force for hydrophobic 
assembly, and act to reinforce the influence of the free energy surface on the association of 
more hydrophilic bodies. 



I. INTRODUCTION 

The hydrophobic effect plays a fundamen- 
tal role in biology and nanoscale engineering. 
Great strides have been made in our under- 
standing of the mechanism of hydrophobic as- 
sembly by elucidating the behavior of solvent 
near interfaces of variable length scalePEH i n 
addition, the hydrophobic effect is crucial to 
understanding the nature of surface friction in 
nanoscale fluid transport G2H2H Surface friction 
characterizes the boundary condition in the con- 
tinuum description for fluid flow across inter- 
faces of varying hydrophobicity. 

Despite recent efforts to improve our un- 
derstanding of the hydrophobic effect, to our 
knowledge the role of hydrodynamic interac- 
tions in hydrophobic assembly has not been con- 
sidered. Within the framework of Brownian dy- 
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namics, hydrodynamic interactions are incorpo- 
rated in the frictional force by means of a resis- 
tance (friction) tensor that is dependent upon 
the degrees of freedom of the Brownian bod- 
iesPS Hydrodynamic interactions are typically 
treated in the simulation of colloidal suspen- 
sions within the continuum description of low 
Reynolds number flowPsHU! Recently, polydis- 
perse colloidal simulations intended to mimic 
cytoplasm reported that hydrodynamic interac- 
tions give rise to the slow diffusion of macro- 
molecules observed in crowded cellular environ- 
ments 

For small separations, the continuum descrip- 
tion is dominated by lubrication effects, which 
yields a friction constant that diverges as two 
bodies come into contact. Such frictional forces 
arise from the assumption that there is a fluid 
element present even as the spacing approaches 
zero. Naturally, this description breaks down 
when molecular scale effects become important 
at small length scales. The continuum descrip- 
tion of the Navier-Stokes equation is expected 
be valid down to length scales of approximately 
1 — 2 nm ) 21 l 24 l Therefore, molecular-scale effects 
should be considered when only a few solvent 
layers separate the solute. 
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The spatially dependent friction tensor can 
be extracted from explicit molecular dynamics 
simulations where all bath degrees of freedom 
are included. This may be achieved in either 
the Brownian limit or when memory effects are 
includedF^H^H To our knowledge, the friction 
tensor between two bodies has only bee n previ- 
ously computed in Lennard- Jones fluidsP^ 36 * 37 ! 
Presently we consider the spatial dependence of 
the friction constant as function of the distance 
between two non-polar spherical bodies in ex- 
plicit water. The spatial friction, in concert 
with the (non-hydrodynamic) potential of mean 
force along the relative direction, facilitates the 
modeling of two-body assembly in the Brown- 
ian limit and an assessment of the importance 
of hydrodynamic and non-hydrodynamic forces. 

This work aims to elucidate the role that 
hydrodynamics plays in hydrophobic assembly. 
We take as our model solute two fullerenes, ei- 
ther two C60s or two C240s. Unsubstituted 
fullerenes are insoluble in waterpS and thus, to 
our knowledge there are no experimental data 
on their transport in water. Nevertheless they 
are a useful model for theoretically probing the 
transport of nanoscale bodies. The hydropho- 
bic collapse of two bodies in explicit water^EH 
and pathways for self-assembly of hydroph obic 
spheres in a coarse grained solvent modePE^ 
have been previously studied. The potential 
of mean force of two C 60 m olecules in water 
has also been computed! 39 * 40 * However, the in- 
terplay between the free energy surface and hy- 
drodynamic interactions has not yet been con- 
sidered. As we will show, solvent fluctuations 
and drying phenomena manifest themselves as 
hydrodynamic interactions within the Brown- 
ian framework. We find that because slow re- 
laxation times are associated with the solvent 
at the drying transition the separation of solute 
and solvent time scales required by a Brownian 
description of hydrophobic assembly may not 
always be satisfied. 

The solute-solvent interactions are described 
by both attractive and purely repulsive po- 
tentials. For the purely repulsive model of 
the fullerene we observe a drying transition^ 
one which is more pronounced for C240 than 



C60, a finding which is in agreement with the 
known length-scale dependence of drying phc- 
nomenaPESl The attractive spheres do not dry 
as they approach each other and water is ex- 
pelled from the inter-solute region by means of 
steric repulsion. We find that the spatial de- 
pendence of the friction constant is quite differ- 
ent in each of these cases. In particular, drying 
phenomena and solvent layering engendered by 
attractive potentials are shown to possess dis- 
tinct signatures in their respective friction ten- 
sor. Hydrodynamic interactions at the molecu- 
lar scale are related to the variability of solvent 
density, fluctuations, and relaxation times in 
the inter-solute region as two bodies approach. 

The rate constant for asse mbly i s computed 
via a Smoluchowski Analysis PE3EI] We find 
that for two body assembly, the inclusion of 
hydrodynamic interactions decreases the rate 
of reaction in the diffusion controlled limit by 
approximately 30 — 40%. In the case of ideal 
hydrophobic assembly (that is when no solute- 
solvent attraction is present) the frictional and 
mean force can be said to have opposite effects. 
When solute-solvent attraction is included, they 
tend to reinforce each ether's impact on the rate 
of assembly. 

This paper is organized as follows. Sec. [TT] 
reviews the formalism of hydrodynamic interac- 
tions and Brownian motion. Simulation details 
are given in Sec. |III| Results are given for the 
friction coefficient on a single body in Sec. |IV| 
and are compared with the Stokes-Einstein rela- 
tion with periodic boundary corrections. In Sec. 
|V| an analysis of the assembly of two non-polar 
solutes within the framework of Brownian mo- 
tion is presented. Discussions and conclusions 
are given in Sec. |VI| 



II. LANGEVIN DYNAMICS AND 
HYDRODYNAMIC INTERACTIONS 

In this work we seek to characterize solvent- 
mediated effects on two spherical bodies. The 
Brownian limit is considered such that the bod- 
ies evolve on a much slower time scale than the 
solvent bath. The dynamics of such a system 
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can be described by an equation of motion that 
includes contributions from three forces, a fric- 
tional force, a mean force, and an uncorrelated, 
Gaussian random force. Due to the symmetry of 
the problem, the relative distance r = \?2 — ?\\ 
between the spheres is the only coordinate nec- 
essary to characterize the approach of the two 
solute bodies. In the present study we concen- 
trate on this direction, and will consider the 
friction along other degrees of freedom in fu- 
ture work. When hydrodynamic interactions 
are included in the frictional force, the Langevin 
equation for the motion of the relative coordi- 
nate is given byj^ 

fir = - VW(r) - ((r)r + R{r, t) , (1) 

where ji is the reduced mass and W(r) is the 
potential of mean force. In the Brownian limit 
the solvent (hydrodynamic) time scale is sepa- 
rable from the motion of the solute and so the 
friction tensor may be computed at fixed so- 
lute positions.^ This is analogous to the Born- 
Oppenheimer approximation, which assumes a 
similar separation of time scales between nu- 
clear and electronic degrees of freedom. The 
(non-hydrodynamic) free energy profile contains 
contributions from both direct solute-solute in- 
teractions and induced solvent- mediated forces. 
In the limit of large friction the term on the 
left-hand side of Eqn. Illmay be neglected. The 
friction coefficient £(r) can be determined from 
the full two-body tensor as reviewed in the Ap- 
pendix. The random force, R(r,t), is Gaussian, 
uncorrelated noise with zero mean and a vari- 
ance that is related to the friction coefficient via 
the fluctuation-dissipation theorem, 

(R(r,t)R(r,0)) = 2k b T((r)6(t) . (2) 

As indicated by Eqn. [l] when hydrodynamic 
interactions are included the friction coeffi- 
cient depends on the distance between spheres. 
If hydrodynamic interactions are ignored then 
C{r) — > Co, where Co is taken to be the friction 
coefficient on the relative coordinate at infinite 
separation where no spatial dependence is as- 
sumed. 

Hydrodynamic interactions that yield the 
spatial dependence of the friction coefficient 



are typically computed by means of approxi- 
mate solutions of the linearized Navier-Stokcs 
equation valid for incompressible flow at low 
Reynolds number. A ubiquitous formulation 
of hydrodynamics frequently utilized in col- 
loidal simulations is Stokesian dynamics! 26 * 27 * 
Stokesian dynamics interpolates between an ex- 
act two-body solution for short range interac- 
tions^ and the long range many-body result 
derived from the Rotne-Prager tensor as ap- 
plied to systems with periodic boundary condi- 
tions.^. Other techniques may also be utilized 
to c ompute co ntinuum hydrodynamic interac- 
tion a 28 * 29 ! 49 * 50 ! and the continuum formulation 
has been recently employed as a framework for 
coarse-grained solvent models.^ In this work 
we will compare our molecular scale results to 
both the expression of Ref. 47 and the Oseen 
tenso r com puted with periodic boundary condi- 
tions E2IH The Oseen tensor is the leading order 
term of the Rotne-Prager expression^. 

The continuum approach breaks down for 
small particle separations, and in this work we 
seek to determine hydrodynamic interactions at 
the molecular scale. The solvent bath is treated 
by explicit Newtonian (microcanonical) molec- 
ular dynamics simulation. All bath degrees 
of freedom and microscopic details are present 
in our simulation. The friction coefficient is 
related to the microscopic fluctuations in the 
linear response regime by means of a Green- 
Kubo relation. This equates the friction ten- 
sor to the correlation function of the fluctua- 
tions of the total force on each Brownian body, 
SFi = Fi-(Fi)m 

OO 

C =T^ f dt J im (^W^(O)) (3) 

=i] Kfjl J N— >oo \ / 


Where £ij is a component of the friction ten- 
sor (see Appendix). Unfortunately, this for- 
mula may only be directly applied in both the 
Brownian limit (M so i ute — > 00) and in the ther- 
modynamic limit where the numb er of so lvent 
molecules, N, approaches infinity.^ 5 * 36 * 54 *. Al- 
though the former requirement may be satisfied 
by fixing the positions of the Brownian parti- 
cles, the latter requirement gives rise to sub- 
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tie complications in a finite simulation. Utiliz- 
ing the techniques developed by Bocquet et Ab. 
Eqn. [3] may be employed to indirectly compute 
the two-body friction coefficient, that is the MD 
estimate of £ can be related to linear combina- 

tions of the submatricies of the friction tensorPS 
Similar complications are present in the compu- 
tation of the friction on a single bodyP^ Further 
details are provided in the Appendix. 



III. SYSTEM SETUP AND SIMULATION 
DETAILS 

In order to gauge the impact of size effects, 
we utilize two fullerenes, C60 and C240, as our 
model non-polar bodies. The two molecules and 
their relative sizes are depicted in panel (b) of 
Fig. [T] The solute-solvent interactions are de- 
scribed by two forms, a Lennard- Jones (LJ) in- 
teraction with parameters <rgs = 0.332 nm and 
e ss = 0.423 kJ/mol, and a Weeks-Chandler- 
Andersen (WCA) truncation of this potential. 
As will be shown below, there is vastly differ- 
ent behavior when the potential is attractive or 
purely repulsive. For the purpose of the friction 
analysis the fullerenes are treated as spheres and 
the friction coefficient is computed on their re- 
spective centers of mass. Based on the choice 
of repulsive core, the estimated Van der Waals 
diameters for C60 and C240 are a = 1.0 nm and 
a = 1.7 nm, respectively. 

The friction is sensitive to periodic bound- 
ary conditions (see Sec. IV), and we utilize as 
large a system as is computationally feasible. 
For runs containing C60 the base box length 
prior to equilibration is 5.0 nm and for C240 it 
is 6.0 nm. In systems containing two fullerenes, 
a set of calculations was run with the bodies 
frozen and set at a fixed distance apart from 
each other. The box length is extended by 
the separation between solute centers of mass 
in the axial direction of the starting configura- 
tion of each computation. The number of water 
molecules utilized to represent the solvent varies 
between approximately 5000— 11000 depending 
on the size of the fullerene and the size of the 
separation. Water is modeled with the TIP4P 



potential^. Snapshots^ of two C240 molecules 
fixed at three different separations are shown 
in Fig. [2] Each system was equilibrated under 
NPT conditions for 2 ns and NVT conditions 
for 1 ns. The barostat of BerendserPS and a 
stochastic velocity rescaling thermostaP^l were 
employed during equilibration runs to maintain 
a temperature of 300 K and a pressure of 1 bar. 
In order to compute the friction coefficient, data 
was collected from 10 — 18 NVE runs of 4 ns 
in length at each fixed distance. All calcula- 
tions were performed using GROMACS version 
4.5.3P3 As the friction is a probe of the solvent 
momentum relaxation in the presence of the so- 
lute bodies, its computation requires strict en- 
ergy conservation in the microcanonical ensem- 
ble. NVE runs were performed in double preci- 
sion. 

The computation of the spatially dependent 
friction coefficient requires constant energy sim- 
ulations and three dimensional periodic bound- 
ary conditions are necessary for comparison 
with long range periodic continuum hydrody- 
namics. Therefore, behavior at the drying tran- 
sition must be studied at constant volume and 
in the absence of an interface. The number of 
solvent molecules we choose is large enough to 
accommodate the study of a dewetting transi- 
tion and although the behavior is modestly per- 
turbed by this choice of ensemble, the qualita- 
tive features of the friction at the critical dis- 
tance for dewetting are not expected to suf- 
fer. As the two solutes are fixed throughout 
all computations, direct solute-solute interac- 
tions are not a component of the simulations 
described above. Direct interactions between 
fullerenes are included in the computation of the 
rate constant in Sec. IV CI and the solute-solute 
potential has the parameters a — 0.35 nm and 
e = 0.276 kJ/mol. 



IV. FRICTION ON A FULLERENE IN 
SOLVENT 

The strength of non-polar attractions play 
an important role in the degree of particle hy- 
drophobicity (see e.g. Ref. |5])- This is readily 
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FIG. 1. (a) The pair-correlation function from the center of C240 (top panel) or C60 (bottom panel) to 
the oxygen site on the water molecule is depicted for the cases when the solute-solvent potential is purely 
repulsive (dashed black curve) or includes attraction (solid red curve), (b) The two fullerenes utilized to 
model non-polar spherical bodies, C60 (left) and C240 (right) are visualized with VMDP The estimated 
diameter for each species is 1.0 nm and 1.7 nm, respectively, (c) Two spherical bodies of diameter a are 
placed a distance r from each center. A cylindrical probe volume (blue lines) is chosen such that its length 
is r — a and its diameter determined by the angle <f> = 90°. This volume can be further partitioned into 
an inner-tube (red lines) whose diameter is determined by the angle 6 = 60° and an outer shell that is the 
difference between the two probe volumes. 



seen in panel (a) of Fig. [Tj where the pair cor- 
relation function of the solvent with respect to 
the center of mass of the fullerene is plotted. 
Whereas the attractive fullerenes exhibit sig- 
nificant structuring of the first solvation shell, 
the water density is depleted near the surface 
of the repulsive bodies. The pair correlation 
function for C60 is in good agreement with pre- 
vious studies of hydration around fullerenes.^ 
Furthermore, depletion becomes more promi- 
nent with increasing solute size, a finding that is 
in agreement with prior theoretical predictions 
and simulation.^E3] 

The friction coefficient of a single C60 and 
C240 molecule is computed utilizing both the 
LJ and WCA solute-solvent potentials. These 
values may be related to the Stokes radius via 
the Stokes-Einstein relation corr ected for peri- 
odic boundary conditions (PBC).^^2l As noted 
in the cited works, PBC have a large impact 
on the measured friction. To leading order this 



correction, as formulated in Ref . 62 , is given by 
the following expression, 



Ciso Cpbc 67T?7i 

where £ = 2.837297, L is the periodic box 
length, and 77 is the shear viscosity. The value 
of the viscosity of the TIP4P model used in this 
work has been taken from the literature.^ 

Once the friction coefficient of the isolated 
sphere, £ iso , has been computed, it can be re- 
lated to the Stokes radius, which is given by, 

Ciso ,_n 

^stokes — 1 [pj 

CUT] 

where c is equal to 4 or 6 for slip or stick 
boundaries, respectively. Studies relating Eqn. 
[5] to the results of molecular dynamics simu- 
lation have appeared in the literature^^U anc j 
subtleties associated with the appropriate hy- 
drodynamic radius and boundary choice should 
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be addressed for a complete understanding of 
how this continuum expression applies at the 
molecular scale. Furthermore, we note that if 
our model solutes were perfect spheres (which 
they are not), then the slip boundary condition 
holds unless solute-solvent attractions are so 
large that solvent molecules become entrenched 
on the surface.^ Presently, we simply assume 
either a stick or slip boundary condition and 
compute the stokes radii of the fullerenes with 
the LJ and WCA solute-solvent potentials. The 
results are presented in Table [I] and compared 
with the estimated van der Waals radii. One 
can see that the stokes radii of the purely re- 
pulsive spheres with slip boundary conditions 
are in very good correspondence with the van 
der Waals radii. As the molecules become more 
attractive, the friction coefficient increases as 
does in this sense the "stickiness" of the sur- 
face. This is in agreement with studies of water 
at interfaces, which have shown that there is 
a correlation between hydrophobicity and the 
friction tangential to the surface E2122I Further- 
more, studies of particle diffusion in a Lcnnard- 
Jones fluid have exhibited similar trend s as a 
function of solute-solvent attraction! 65 ! 66 ! 

In the following analysis, we determine the 
friction coefficient along the relative coordinate 
in units of Co = Cpbc/2 for each respective 
fullerene and solute-solvent potential. The two 
body simulations utilize a periodic box of di- 
mensions (L + d) x L x L and therefore this 
choice accounts for the impact of periodic con- 
ditions. Changes relative to Co are a result of 
hydrodynamic interactions and not PBC arti- 
facts. 



V. TWO-BODY RELATIVE FRICTION 

The friction coefficient on the relative coor- 
dinate has been computed as a function of the 
distance between two spherical bodies where the 
solute-solvent interaction is either attractive or 
purely repulsive. We find starkly different be- 
havior in each will be discussed below. 
In both cases, there is significant deviation from 
continuum hydrodynamic predictions at small 



separations. In Fig. [2] we plot the relative 
friction coefficient as a function of the sphere 
separation r — a for C240 where solute-solvent 
interactions are purely repulsive or include at- 
traction. Two expressions from continuum hy- 
drodynamics, the two-body result valid at small 
separations and the result of the Oseen tensor 
corrected for periodic boundary condi tions that 
is valid at long range, are also plotted. ^ 7 * 50 * 52 * 53 ! 
One can see that although the friction coeffi- 
cient for r — a > 1 nm, approaches the long- 
range continuum result, the short-range value 
does not diverge at the point of contact as pre- 
dicted by the expression of Ref. |47] Instead, 
the frictional profile exhibits distinctive features 
of molecular-scale effects. For purely repulsive 
bodies, the friction coefficient peaks at a dis- 
tance that correspond s to the critical distance 
for drying (see Sec. 



VA) 



This distance in- 
creases with solute length scaled and such dry- 
ing phenomena is not encompassed by contin- 
uum hydrodynamics. In the case of bodies that 
strongly attract solvent, drying does not occur, 
and water is expelled by steric repulsion. This 
result is dependent upon the chosen value of 
655. More weakly attractive spheres can ex- 
hibit dewetting. As discussed in Sec. |VB| dis- 
tinctive molecular-scale behavior in which the 
friction coefficient exhibits signatures of solvent 
layering is observed as solvent is expelled from 
the inter-solute region. Furthermore, we note 
that the value of the long-range molecular scale 
results plotted in Fig. [2] appear to approach 
the description given by the periodic Oseen ten- 
sor. However, data points at larger separa- 
tions would be necessary to make more defini- 
tive comparisons. 

The occupancy and fluctuations of solvent 
molecules in a particular region of space may be 
monitored in a probe volume. This technique is 
a standard tool to study the hydrophobic effect 
(see e.g. Ref. |2J. Nested cylindrical probe vol- 
umes are presently employed to study the water 
density and fluctuations in the inter-solute re- 
gion. The probe cylinder is partitioned into an 
outer shell and an inner-tube. Due to the cur- 
vature of the hydrophobic surfaces, the density 
may be significantly different in the two shells, a 
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Stokes radius (slip) 


Stokes radius (stick) 


Estimated vdW radius 


C60 


WCA 


0.53 inn 


0.35 nm 


0.50 nm 


C60 


LJ 


0.66 nm 


0.44 nm 


0.50 nm 


C240 


WCA 


0.84 nm 


0.56 nm 


0.85 nm 


C240 


LJ 


1.06 nm 


0.70 nm 


0.85 nm 



TABLE I. Table of computed Stokes radii for C60 and C240 with differing solute-solvent interactions. 



finding that is particularly true in the case of the 
fullerenes which attract solvent. A diagram of 
the boundaries of the probe volume is sketched 
out in panel (c) of Fig. [lj The length of the 
cylinder is determined by the separation of the 
surfaces of bodies. The diameter of the cylin- 
der is formed by a line tangent to the effective 
spherical surface of the fullerene and opposite 
an angle of 90°. The inner-tube is bounded in 
the same fashion except the angle opposite the 
tangent segment is 60°. 



A. Purely repulsive solute-solvent interactions 

In this section the results when only repul- 
sive interactions with the solvent are present 
is discussed. The frictional profiles for C60 
and C240 are given alongside their respective 
solvent-induced potentials of mean force (PMF) 
in panels (a) and (b) of Fig. [3j It can be seen 
that the hydrophobic collapse along the relative 
coordinate is essentially barrierless. This find- 
ing is in agreement with previous simulations 
of the assembly of two hydrophobic bodies in 
bo th coars e-grained and explicit solvent mod- 
els jMHEH The density and fluctuations of water 
in the inter-solute regions are depicted in Fig. |4j 
The onset of drying is at larger distances with 
increasing solute size, an observation which is in 
agreement with prior work.™ The drying tran- 
sition for the approach of two C240 molecules 
occurs when approximately two layers of wa- 
ter would be present in the inter-solute region 
if only steric effects were considered. Upon in- 
spection of Fig. [4j it is apparent that the ratio 
of density of solvent to its bulk value is less than 
unity even at large separations when the region 



is wet, an indication that water is depleted near 
the surface of the hydrophobe, a result in agree- 
ment with the pair-correlation function shown 
in Sec. El 

The fluctuations of water in the probe volume 
are plotted in the lower panels of Fig. [4] The 
chosen metric of fluctuations / (N) for a 

bulk fluid is proportional to the isothermal com- 
pressibility in the limit of large probe volumes, 
which naturally does not presently apply. The 
fluctuations peak as the drying transition oc- 
curs. These findings are consistent with the un- 
derstanding of dewetting phenomena that have 
been put forth in previous studies.™ 

In order to gain an understanding of how 
molecular-scale effects manifest themselves in 
the spatially dependent friction coefficient, we 
now show how the value of the friction coeffi- 
cient relates to drying phenomena. The fric- 
tional profile is depicted in panel (a) of Fig. [3] 
and is strongly peaked at the dewetting tran- 
sition. For larger solutes such as C240, this 
may be defined as the separation where the 
inter-solute region fluctuates between "wet" and 
"dry" states.™ The histo gram of the number of 
waters in the probe volume at r c is shown in 
the top panel of Fig. [5j There is a bimodal dis- 
tribution where the difference between maxima 
is ss 10 water molecules. The distributions at 
r = r c ± 0.1 nm are also shown, and at these 
separations the inter-solute region is predomi- 
nantly wet or dry. Although drying occurs as 
two C60 fullerenes approach each other, a bi- 
modal distribution as exhibited in Fig. [5] is less 
apparent. In addition, the critical separation 
is near the inflection point of the density curve 
and the maximum of the solvent fluctuations as 
determined by the chosen metric. The point at 
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FIG. 2. The friction coefficient is plotted as a function of the distance between two C240 molecules when 
the WCA potential and LJ potential is utilized to model solvent-solute interactions (black dotted line with 
circles and red dashed line with squares, respectively). The two-body (continuum) hydrodynamics result of 
Ref. [47] (solid blue line) and the periodic Oseen tensor (dashed green line) are plotted against the molecular- 
scale computation. The three right-hand panels depict snapshots at separations (increasing from top to 
bottom) where, in the case of the WCA spheres, the inter-solute region is dry and the friction coefficient 
is low (violet circle), at the dewetting transition where the friction peaks coefficient (orange circle), and 
for large distances where the friction coefficient decays towards the baseline result (green circle). Water 
molecules within 0.6 nm of the fullerene surface are depicted by a ball-and-stick representation. All others 
are represented as lines. The snapshots were visualized with VMDP 



which the friction coefficient peaks is marked in 
the graphs plotted in Fig. [4] 

In addition to large static fluctuations (see 
Fig. [2]), slow relaxation times are also in- 
dicative of the drying transition. The au- 
tocorrelation function of the fluctuations in 
the number of water molecules in the probe 
volume, (SN(t)6N(0)), may be decomposed into 
a product of the variance (SN 2 ) and some time 
dependent function h(t) that characterizes the 
relaxation, 

(SN(t)8N(0)} = (8N 2 ) h(t) . (6) 

The normalized time correlation function h(t) 
and the variance are plotted in the bottom panel 
of Fig. [5j One can see that, in addition to 



a large variance, slow relaxation times are as- 
sociated with collective water dynamics at the 
critical distance for drying. In this way both 
static and dynamic terms contribute to the sol- 
vent behavior that engenders a large friction co- 
efficient. The relaxation of the correlation func- 
tion is markedly faster when r = r c ± 0.1 nm, 
and the friction coefficient at these separations 
is lower by a factor of « 4. The computation 
of the hydrodynamic interactions is a probe of 
the momentum relaxation of the solvent (see 
Appendix) and slow collective motion of water 
molecules exiting and entering the inter-solute 
region gives rise to, in part, the large friction 
that is observed. However, as the time scale 
of solvent relaxation approaches that of the so- 
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lute, the separation of time scales essential to 
the (Markovian) Langevin description will not 
hold. Such possibilities will be discussed in Sec. 

ED 

The distribution of water molecules in the 
same size probe volume used at r = r c is given 
near the surface of a single C240 and in the 
bulk in Fig. [6] The depletion of water density 
at the surface is indicated by the differences in 
the bulk and single surface distributions. Fur- 
thermore, the enhancement in fluctuations near 
the surface is indicated by the broader distri- 
bution. These findin gs are in broad agreement 
with prior work.LLIUllThe faster relaxation of the 
number correlation function is seen in the bulk 
and near the surface. This indicates not only 
the slow collective motions associated with the 
bi-stability of the drying transition, but also the 
more general effect of confinement between the 
two solutes lengthens the relaxation times in the 
probe volume. 



B. Attractive solute-solvent interactions 

Given the high density of carbon sites and 
the value of ess chosen for the Lennard- Jones 
potential, the solutes possess a sizable affinity 
for water when attractive interactions are in- 
cluded. The solvent density is plotted in Fig. 
[7[ The curvature of the surface has a signifi- 
cant impact on water density and fluctuations. 
For small separations where solvent in the inter- 
solutc region can experience attraction from two 
bodies, the density is enhanced in areas in the 
outer probe region where steric repulsion does 
not dominate. High densities are accompanied 
by large fluctuations of the number of water 
molecules in the outer shell. Due to surface 
curvature, water is first squeezed out from the 
inner probe volume as the density decreases in 
this region as the second layer is vacated, and 
subsequently, the density falls to zero as the fi- 
nal layer of water is removed. In analyzing the 
trends in the friction coefficient and potential 
of mean force, the inner probe volume is more 
indicative of solvent layering. 

The induced potential of mean force depicted 



in the panel (d) of Fig. [3] exhibits minima at 
separations where one and two water layers are 
stable at the narrowest portions of the inter- 
solute region, as indicated by monitoring the 
inner probe volume. As the spheres come into 
contact there is a maximum in the solvent in- 
duced PMF. Prior calculations of the PMF for 
a pair of C60 molecules are in qualitative agree- 
ment.^! Prior computations of the PMF as two 
graphene plates approach exhibit barriers for 
the removal of each solvent layerjS3 although 
the curvature of the spheres gives rise to differ- 
ent behavior than plates at small separations. 
Namely, as the spheres are curved, the inter- 
solute region never fully desolvates and the in- 
duced mean force opposes assembly until the 
point of contact. 

Upon examination of the spatial dependence 
of the friction coefficient (panel (c) of Fig. [3]), 
evidence of layering is exhibited by the non- 
monotonic behavior of the friction coefficient for 
small separations r — a < 0.5 nm. The fric- 
tion coefficient peaks as the PMF rises and the 
second layer of solvent is removed. At smaller 
separations, the friction coefficient falls at the 
minimum of both the induced PMF and sol- 
vent fluctuations (r — a = 0.3 nm), and then 
as the final water layer vacates the inter-solute 
region both the PMF and friction coefficient 
rise sharply, with the size of the latter abat- 
ing rapidly as solvent density in the inner re- 
gion decreases. These trends are also witnessed 
in the water statistics. The point at which the 
friction coefficient is maximum corresponds to 
a point where the density is near a maximum, 
and the fluctuations are large and increasing. 
The values of the density and fluctuation at the 
separation of maximum friction coefficient are 
marked on the curves presented in Fig. [7] In the 
case of C240, these trends are more clearly ob- 
served upon study of the smaller probe volume. 
It stands to reason that in the case of attrac- 
tive spheres, the friction exhibits some direct 
relation to both the density and fluctuations of 
solvent in the inter-solute region. 

The spatial dependence of the friction coef- 
ficient is in stark contrast to the predictions 
of continuum hydrodynamics in the short-range 
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FIG. 3. The spatial friction coefficient (panel a) and solvent-induced potential of mean force (panel b) as 
two C60 (red curve) and two C240 (black curve) approach each other when the solute-solvent interactions 
are purely repulsive. As expected, there is a greater driving force for assembly with increasing solute 
length-scale. Panels (c) and (d) depict, respectively, the spatial friction and solvent-induced potential of 
mean force as two C60 (red curve) and two C240 (black curve) approach each other when the solute-solvent 
interactions include attraction. In this case, no drying is observed and there is a barrier for desolvation of 
two water layers as the two bodies approach each other. The friction coefficient is not given at the point of 
contact for the WCA spheres due to the substantial uncertainty in the correlations of the force fluctuations 
when the mean force is large. 



limit, where the friction coefficient is predicted 
to monotonically increase and diverge as the 
bodies are brought into contact. This discrep- 
ancy is engendered by the layering and granu- 
larity of liquids over small length scales. The in- 
crease in friction coefficient may have its molec- 
ular origins in the confined rattling in the inter- 
solute region as the water is squeezed out as is 
indicated by the large fluctuations in the probe 
volume (Fig. [7|. It is also observed that the 
spatial friction coefficient plotted in panel (c) of 
Fig. [3]scales with particle size, that is the curves 
for C60 and C240 nearly superimpose upon each 
other when the difference in effective molecular 
radius is taken into account. 



C. Diffusion controlled rate 

In the high friction limit, a radial Smolu- 
chowski equation describes the motion along the 
relative coordinate, r. The diffusion-controlled 
rate of reaction can b e computed from the fol- 
lowing expressio n 1 33 ! 41 ! ^: 

(7 

where a is the contact diameter of the Brownian 
bodies and, 



(8) 
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FIG. 4. The relative density of water (top row) and the ratio of the variance to the average of the number 
of water molecules (bottom row) in the probe volume is plotted for C240 (right column) and C60 (left 
column) when the solute-solvent interaction is purely repulsive. The probe region is divided into an outer 
shell (red dashed line) and an inter-tube (blue dashed line). The squares indicate the separation at which 
the friction coefficient is maximum. 



is the ratio of the friction coefficient to the pair- 
correlation function g(r) — e~ w ( r ^ kbT , where 
W(r) is the total potential of mean force that 
includes contributions from both direct and (in- 
duced) solvent-mediated interactions. The full 
PMF is plotted in panels (b) and (d) of Fig. [8] 
and for C60 is in good agreement with prior 
workP^l The function Q(r) is plotted in pan- 
els (a) and (c) of Fig. 8] for all combina- 
tions of fullerene sizes and interaction poten- 
tials presently considered. It is compared with 
its value if hydrodynamic interactions are ne- 
glected, that is when £(r) = Co, where Co is one 
half the friction coefficient on the single body 
under periodic conditions (see Sec. IV I. 



In the case of purely repulsive solute-solvent 
interaction, Q(r) has a markedly different form 
that depends on whether or not hydrodynamic 
interactions are included. In particular, due to 
the fact that the PMF is virtually barrierless 
when solute-solvent interactions are purely re- 
pulsive, hydrodynamic interactions give rise to 



a maximum in Q(r), as exhibited in panel (a) 
of Fig. [5J In the case of bodies that attract sol- 
vent, Q(r) is also larger when hydrodynamic in- 
teractions are included. However in contrast to 
ideal hydrophobes, the inclusion of spatially de- 
pendent friction coefficient tends to "reinforce" 
the extrema engendered by the PMF. Consis- 
tent with the frictional profiles presented in Fig. 
[3j the impact of hydrodynamic interactions is 
greatest at separations where drying or desol- 
vation occurs. The function Q(r) is found to 
increase at small separations, r — a < 0.1 nm, 
as direct repulsive interactions dominate in this 
region. 

In order to evaluate Eqn. [7J the spatial fric- 
tion coefficient is approximated for large separa- 
tions via a switching function 68 that smoothly 
reduces its value to a plateau of Co at a dis- 
tance 0.5 nm further than the farthest separa- 
tion for which the friction coefficient was com- 
puted. Although the computed rate constant 
exhibits weak dependence on the choice of cut- 
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FIG. 5. (Top panel) The distribution of number of 
water molecules present in the probe volume at the 
critical distance for drying, r c — a = 0.6 nm (black 
curve) and at r — r c ±A where A = 0.1 nm (red and 
blue dashed lines). (Bottom panel) At the critical 
distance there is a slow relaxation observed in the 
autocorrelation function of the number of molecules 
in the probe volume as compared to states that are 
primarily dry or wet. The inset depicts the variance 
of the number of waters in the probe volume as a 
function of A. All data in this figure pertains to 
the C240 system when solute-solvent interactions 
are purely repulsive. 



off, this does not alter the observed trends. We 
note that Co depends on the periodic bound- 
ary conditions adopted (see Sec. IV), and, 



moreover, due to the small relative separations 
studied in our work, our calculations do not 
capture the long-range decay of the hydrody- 
namic interactions. Thus the rate constant for 
association presently computed accurately re- 
flects the short-range hydrodynamic interaction 
at the molecular scale but not the contribution 
of the long range hydrodynamic interactions. 

If both the potential of mean force and hydro- 
dynamic interactions are neglected the diffusion 
controlled rate constant reduces to the familiar 




FIG. 6. The distribution of water molecules in a 
probe volume of length r c — a near the surface of 
a single C240 molecule (solid black curve) and in 
the bulk (dashed red curve) when the solute-solvent 
potential is purely repulsive. The inset depicts the 
autocorrelation function of the number of waters in 
the probe volume in the respective environments. 



expression, 



k 



Airk^Ta 



(9) 



The ratio of the rates to ko computed by means 
of Eqn. [7] with and without the inclusion of 
hydrodynamic interactions or the potential of 
mean force are given in Table [TTJ The values 
of ko are nearly independent of molecular di- 
ameter, which is a consequence of the Stokes- 
Einstein relation and the direct proportionality 
of friction coefficient to solute diameter. It can 
be seen that if the potential of mean force is 
considered and hydrodynamic interactions are 
not included (that is £(r) = Co), then the rate 
of "ideal" hydrophobic assembly (when the so- 
lute only excludes volume) is increased. In the 
case of the attractive solutes, there are bar- 
riers to assembly and an attractive basin at 
short range that serve to only modestly impact 
the rate. The inclusion of hydrodynamic in- 
teractions reduces the diffusion controlled rate 
by approximately 30 — 40% from its respective 
value across all systems presently considered. 
If the potential of mean force (both direct and 
solvent-induced) is not included in the compu- 
tation of the rate constant, then one can clearly 



13 




FIG. 7. The relative density of water (top row) and the ratio of the variance to the average of the number 
of water molecules (bottom row) in the probe volume is plotted for C240 (left column) and C60 (right 
column) when the solute-solvent interaction includes attraction. The probe region is divided into an outer 
shell (red dashed line) and an inner-tube (blue dashed line). Triangles and squares denote the separation 
at which the friction coefficient is maximum. 



see that hydrodynamic interactions reduce the 
rate of association almost uniformly across the 
four systems studied, this despite the different 
physics in their separate behaviors (see Fig. [8]) . 
This is probably due to the fact that the main 
quantitative impact of molecular-scale hydrody- 
namics is to increase the value of Q(r) at sep- 
arations slightly larger than required for desol- 
vation of the inter-solute region. 

Continuum hydrodynamic theory has been 
utilized to compute the rate of assembly of 
two bodiesPHESESI As is well known, the 
strong divergence of the friction coefficient ex- 
hibited by the short-range lubrication expres- 
sion utilized in Stokesian dynamics and plot- 
ted in Fig. [2] prohibits contact and thus the 
comp uted r ate constant from Eqn. [7] would be 
zer0 [26127147J rpj^ lubrication limit for the case of 
slip boundary conditions as derived by Wolynes 
and DeutclpS is integrable and hydrodynamic 
interactions were found to decrease the rate con- 
stant by ~ 30%. Recently, hydrodynamics were 



included by means of the Rotne-Prager tensor 
in the evaluation of protein-protein association 
rates and found to decrease the rate of associ- 
ation by w 35 - 80%P Although these results 
generally appear in agreement with our work, 
as noted above, the present estimates only ade- 
quately gauge the impact of short-range hydro- 
dynamic interactions and therefore do not read- 
ily lend themselves to detailed comparisons. 

As stated above, if the spatial dependence 
of the friction coefficient is neglected, the rate 
constant for reaction is larger than ko for ideal 
hydrophobes, as the solvent induced potential 
of mean force is attractive. This rate constant 
is reduced by hydrodynamic interactions which, 
in some sense, capture the timescale associated 
with dcwctting fluctuations (see Sec. VA). As 
drying fluctuations are relatively slow (see Fig. 
[5]) this raises the question of whether or not 
the separation of timescales within the Brow- 
nian limit is good description of systems near a 
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FIG. 8. The ratio of the friction coefficient to pair-correlation function, Q(r) is plotted with and without 
the inclusion of hydrodynamic interactions (HI) as two C240 (black curves that are dotted with squares 
and dashed, respectively) and two C60 (red curves that are dotted with circles and dashed, respectively) 
molecules approach when the solute-solvent interactions include attraction (panel a) and for the case when 
the solute-solvent interaction are purely repulsive (panel c). Q(oo) = £o is defined as the value of Q for large 
separations when the potential of mean force and hydrodynamic interactions are negligible. The potential 
of mean force including the direct interactions between a pair of C240 and C60 molecules with solute-solvent 
interactions that are purely repulsive (panel b) and include attraction (panel d) is also depicted. 



drying transition. 



VI. DISCUSSION AND CONCLUSIONS 

The spatial dependence of the friction coef- 
ficient in the Brownian limit has been com- 
puted as two non-polar bodies come into con- 
tact. The system is analyzed within the frame- 
work of Brownian dynamics and hydrodynamic 
interactions are included via explicit molecular 
dynamics simulation. We find the friction coef- 
ficient deviates from continuum hydrodynamic 
predictions at small separations and dramati- 
cally depends on the nature of solute-solvent 
interactions. For purely repulsive spheres, we 
find that the friction coefficient peaks at the 
critical distance for dewetting and decreases as 
the inter-solute region dries. For attractive so- 



lutes water is expelled by steric repulsion and 
the effects of solvent layering are apparent in the 
non-monotonic dependence of the friction coef- 
ficient on separation. The variation of friction 
coefficient with separation is due solvent occu- 
pation and fluctuations in the inter-solute re- 
gion. Large solvent fluctuations and slow relax- 
ation times are associated with increasing fric- 
tion coefficient, and in this way both static and 
dynamical solvent-mediated effects contribute 
to the frictional force. Under certain circum- 
stance slow solvent fluctuations in the near dry- 
ing regime may give rise to non-Markovian ef- 
fects. 

In their study of hydrophobic assembly, 
Willard and Chandler found that not only the 
relative separation, but also solvent degrees of 
freedom, namely the occupancy of the inter- 
solute region, are necessary to characterize the 
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k/k (HI) 


k/k (No HI) 


k/k (No PMF) 


k ( 10 10 M" 1 s" 1 ) 


C60 


LJ 


0.66 


1.02 


0.68 


1.16 


C240 


LJ 


0.60 


0.95 


0.74 


1.10 


C60 


WCA 


0.93 


1.27 


0.72 


1.54 


C240 


WCA 


0.99 


1.37 


0.64 


1.51 



TABLE II. Table of the ratio of the diffusion controlled rate constants with and without the inclusion of 
hydrodynamic interactions or the potential of mean force to the quantity fco = 4nki,Ta/(o 



process.^ Such behavior has also been observ ed 
in the assembly of hydrophobic plates pSEEH As 
our analysis is in the framework of Brownian 
dynamics, only solute degrees of freedom are ex- 
plicitly treated in the rate calculation presented 
in Sec. |V C[ However, solvent fluctuations en- 
gender the peak in the spatially dependent fric- 
tion coefficient observed at the dewetting tran- 
sition. In this way, drying phenomena manifests 
itself when molecular-scale hydrodynamic inter- 
actions are included in the Brownian limit. 

The Brownian description is valid if the mo- 
tion of the solute occurs on a much slower 
time scale than the solvent dynamics. We find 
that relatively slow solvent motions are present 
at the critical separation for dewetting, which 
raises questions about the validity of the Marko- 
vian assumption. We found the correlation time 
for water fluctuations in the inter-solute region 
at the critical separation to be ps 25 ps for 
C240. (see Fig. [H). Prior work has shown that 
the timescale for hydrophobic assembly depends 
on not only the size of the s olutes but also on 
the nature of the interactions f^EH and the time 
scale for solvent fluctuations can be much longer 
than 25 ps. In the case where solvent fluctua- 
tions cannot be considered fast with respect to 
solute diffusion, one may consider alternative 
formulations that include dynamical disorder^, 
and the use of techniques for extracting the 
spatially dependent memory kernel within the 
framework of the generalized Langevin equa- 
tionPBSSl Naturally, the time scales will depend 
upon the size of the solute, and future studies 
may be designed to more fully probe this obser- 
vation. 

Recently, hydrodynamic interactions have 



been shown to significantly decrease in vivo dif- 
fusion in cellular environments.^ As the present 
computation is limited to two bodies in a sol- 
vent bath, one can only draw limited parallels to 
behavior in crowded systems. With this caveat 
in mind, we find some correspondence in our 
study where we find that inclusion of hydrody- 
namic interactions leads to a decrease of the dif- 
fusion controlled rate constant for assembly by 
approximately 30 — 40%. For barrierless "ideal" 
hydrophobic assembly along the relative coor- 
dinate, hydrodynamic interactions introduce a 
frictional "barrier" that retards the rate of as- 
sembly (see Fig. [8]), whereas in the case of at- 
tractive solvents the main effect is to enhance 
the maxima engendered by the potential of 
mean force. In this way, the interplay between 
hydrodynamic interactions and the free energy 
surface strongly depends on solvent-solute cou- 
pling. Future work will be aimed at exploring 
these insights and applying them to study dif- 
fusive phenomena in nanoscopic systems. 
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Appendix A: Calculation of the spatially 
dependent friction coefficient 

1. Two-body friction tensor 

In this work, we study the spatially depen- 
dent friction coefficient along the relative dis- 
tance between two bodies, r = |r*2 — The 
complete friction coefficient tensor is related to 
the motion of all 3N degrees of freedom. The 
tensor relates the frictional force to the particle 
velocities. For N = 2 this can be expressed as, 



C (r) C (r) 

=11 =12 

C (r) C (r) 

=21 =22 




(Al) 



where 8F is the deviation of the force on the 
sphere from its mean value. The friction coef- 
ficient tensor £(r) has dimension of 6 x 6 and 
may be decomposed into four blocks that cor- 
respond to self and cross interactions between 
the bodies. Each submatrix is a diagonal 3x3 
matrix, 



(A2) 



if the coordinate system is defined such that 
one direction is parallel and two directions are 
perpendicular to the line of centers. The two 
spheres are identical so by symmetry £ (r) = 

C (r) and C (r) = C (r). The friction coeffi- 

=22 =12 =21 

cient along the relative coordinate, r = r 2 ,|| — 
ri ii , may be obtained from manipulation of 
Eqn. | Al | The following relation can then be 
extracted, 




CrelM 



2 (WO 



Ci2,||(r)) 



(A3) 



and associated to the Langevin equation given 
by Eqn. [TJ 
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FIG. 9. The time integral of the self (black curve) 
and cross (red curve) force-force autocorrelation 
functions and the sum of these two functions (green 
curve) is plotted as two C240 molecules approach 
each other and interact with the solvent purely by 
means of repulsion. The top panel, center, and bot- 
tom panel depict this quantity at different separa- 
tions. One can see that the self and cross auto- 
correlation functions plateau at a non-zero value, 
though their sum decays to zero, in agreement with 
the finding of Bocquet et a/P^ The friction coeffi- 
cient is plotted in units of (kj/mol) ■ (ps/nm 2 ). 



2. Review of techniques relating the friction 
coefficient to molecular dynamics simulation 

Here we review the techniques developed in 
Refs. 1551 and I5B1 for the extraction of the Brow- 
nian friction coefficient from molecular dynam- 
ics simulations. First consider a single Brow- 
nian sphere of mass, M, in a bath of N sol- 
vent molecules of mass, m. This comprises an 
isolated system. Next we consider the limit, 
M — > oo, in which the solute particles are fixed. 



17 



In this case, the total solvent momentum is not 
a conserved quantity. As shown below, the com- 
putation of the friction coefficient is essentially 
a probe of the total momentum relaxation. 

We begin by relating the rate of change of 
the total solvent momentum P(t) to the force 
acting on the solute T(t) by means of Newton's 
third law. 



T{t) = -P{t) 



(A4) 



To simplify the notation we only consider one 
dimension, although the expressions for three- 
dimensions are readily obtainable. Next, con- 
sider the integral that in the t — > oo and N — >• oo 
limit yields the Green Kubo relation for the fric- 
tion coefficient. This can be expressed in terms 
of the total solvent momentum: 



i 



dr^(0)J-(-r)) (A5) 



P(t)P(0)) (A6) 



o 

1 



Onsager's principle linearly relates the force act- 
ing on the solute to the total solvent momentum 
in the long time limit: 



T{t) = 



P(t) = 



Nm 
C 



Nm 



P(t) 
P{t) 



(A7) 
(A8) 



Where E qn. |A8| arises from combining Eqn. |A4 
and Eqn. |A7| Utilizing Eqn. |A6| it can readily 
shown that £(t) has the simple form: 



C(t) = (e-< t/Nm . 



(A9) 



Both the t — > oo and N — > oo limits must be 
applied in order calculate the friction coefficient 
in the linear response regime (See Eqn. [3j). If 
the time limit is taken first, as is necessarily 
the case when computing the property in a sim- 
ulation of finite size, then C(oo) — > 0. How- 
ever, if the thermodynamic limit is taken first 
(N —> oo), then a finite and correct value for 
the friction coefficient may be recovered. Prior 



5 
4 

3 
2 

<% 
-1 

-2 

-3 

-4 



■q 

8- • rjnS : : 8 : ■ d '• '• S '■ '•© ■ * - ■■& • <s- ■■ §•■ ■■ 
8 g 11 ? 8." ^ 6 8 8 ° 8 



ii.ll 



^12.1 



5 - 
4 - 
3 - 
2 '- 
j^p 1 J" ' 
- ?. 
»J> .1 - 

-2 - 
-3 - 
-4 - 
-5 - 

it 



0.3 0.6 0.9 1.2 
-> ij ' ' 1 ' ' 1 ■- 



■ Q 

■8 □■□■on ■ '•' -S---0 o-.S-.-.-./.-.vfi.- 

= □ □ □□ g: § ..g :B- 8 

■o o- 
o 

o 



Q ■■■■-{' 



0.6 0.9 1.2 1.5 

r - a (nm) 

FIG. 10. The unique components of the friction 
coefficient tensor as two C240 molecules approach 
each other and interact with the solvent with at- 
tractive (top panel) or purely repulsively (bottom 
panel) forces. The self term parallel (black circles) 
and perpendicular (blue squares) to the line of so- 
lute centers is plotted alongside the cross friction 
coefficient in the parallel (red circles) and perpen- 
dicular (green squares) direction. 
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work has shown that instead of directly apply- 
ing the Green-Kubo relation, the friction coef- 
ficient may be recovered by probing the relax- 
ation of C(t). Presently this was achieved by an 
analysis of the Laplace transform of Eqn. 
developed in Ref. [35] 



In the case of the two body friction coefficient 
we begin by defining the integral of the naive 
MD estimate for each submatrix of the friction 
coefficient tensor that one arrives at from the 
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direct application of the Green-Kubo relation: 

* 

C " )(t) = f^I dT ( SF l® 5Fl W (A10) 
o 

i 

C ^ )(<) = f^rJ dT (^W^(o)) (All) 



The superscript N denotes the number of sol- 
vent molecules present in the simulation. When 
two Brownian bodies are present, the solvent 
momentum relaxation is related to the sum of 
the fluctuations of forces on the two spheres, 
P(t) = -(SFi(t) + 5F 2 (t)). Analogous to case 
of a single body in a bath reviewed above, when 
the long time limit is taken for finite N, it was 
shown in Ref. 1551 that 

oo 

dJ°(oo) + Cg°(oo) = - / dr(p(t)5F 1 (0)) 

o 

= (A12) 



This spurious result is due to finite size effects 
and is demonstrated numerically in Fig. [9j 
Although the measured force-force correlation 
functions of Eqn. |A10| and Eqn. |A11| plateau 
to non-zero values, their sum does indeed de- 
cay to zero. Bocquet et al. proceed to relate 
the integrals of the force-force autocorrelation 
functions with their proper N — > oo limits. De- 
tails are given in Ref. This analysis yields 
relations of the plateaus of the components of 
Q N ' (i) to a linear combination of components 
of the friction coefficient tensor ( that correctly 
approach the thermodynamic limit. 




(A13) 

It can be readily seen that the component of 
the linear combination parallel to axis of centers 
can be identified with the friction coefficient in 
the relative directions that is reported in Fig. 
[3j In order to obtain the full friction coefficient 
tensor, one more relation must be identified, 



2 

NmkbT 




which characterizes the relaxation of the sum of 
the force-force correlation functions plotted in 
Fig. [9] to zero. The resulting friction coefficient 
tensor for the C240 fullerene model with purely 
repulsive and attractive solute-solvent interac- 
tions is plotted in Fig. [TUj It can be seen that 
the friction coefficient is primarily perturbed 
in the direction of the relative coordinate, and 
the perpendicular components of the self sub- 
matrix (Cn,||) have limited spatial dependence 
and their values are closer to the single body 
result. Concurrently, the value of Ci2,_L remains 
near zero. 

There is less than a 10% error in the estima- 
tion of the friction coefficient provided by Ref. 



In 4ln , (A14) 

1 

1361 Estimates of error have been obtained from 
the variance of the measured plateau or from the 
difference between elements of the full friction 
coefficient tensor that are equivalent by symme- 
try. 
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